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Abstract 

We report the results of a search for B~ J/t/jAp based on a data set of 78 fb^^ data collected at 
the T(45') resonance with the Belle detector at the KEKB asymmetric e+e^ collider. No substantial 
signal is found, and we set the branching fraction upper limit B{B~ — J/%l)Kp) < 4.1 x 10~^ at 
90% confidence level. 

PACS numbers: 13.25.Hw,14.40.Gx,14.40.Nd 
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The inclusive decay oi B ^ J /ip + X has been studied by CLEO[l|, Belle [2], and recently 
BaBar|3|. The J/ip momentum spectrum in the e~^e~ center of mass frame is consistent 
with the distribution predicted by non-relativistic QCD calculations except for an excess 
in the low momentum region which has been observed by all of the above experiments. The 
excess below 0.8 GeV/c corresponds to a branching fraction of 4 x 10~^ — 6 x 10~^. 

In order to explain this excess, several theoretical hypotheses have been proposedj^0,0|. 
One of them is that the excess arises from B J/ipKp decaysj^; this possibility can also 
be inferred from the rather large branching fraction of — ^ ApX, (2.3 ± 0.4)% j^. At 
the quark level, B^ J/ipKp can be described as {udud) produced by gluon emission 
from the Cabibbo favoured h ccs diagram. The decay rate could be enhanced by an 
intermediate exotic state such as a bound state of A and p, J / ip and A, or J/ and p. In this 
case, the momentum distributions of the daughter particles will exhibit some characteristic 
enhancements. Thus, searching for B^ —>■ J/ipAp helps to understand the source of the 
excess at low J/ ip momentum and to find intermediate states. The BaBar collaboration has 
recently reported results of a similar search 9]. 

In this paper we report on the study of B^ —>■ J/tpAp. The analysis is based on a data 
sample of 78 fb~^ accumulated at the T(4S') resonance with the Belle detector at the KEKB 
8 GeV and 3.5 GeV asymmetric collider [13]. 

The Belle detector consists of a three-layer silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), an array of aerogel threshold Cerenkov counters (ACC), time- 
of-flight scintillation counters (TOF), a CsI(Tl) crystal electromagnetic calorimeter (ECL), 
a 1.5 T superconducting solenoid coil and an instrumented iron-flux return for muon and 
Ki detection (KLM). The detector is described in detail elsewhere 

In this analysis, we use the decay chain: B~ — ^ J/ipAp, J/ip {l=e,fi), and 

A — > pvr". Inclusion of charge conjugate states is implied throughout this paper. To sup- 
press continuum backgrounds we require the ratio of the second to zeroth Fox- Wolfram 
momentsfH] to be less than 0.5. 

To remove charged tracks that are badly measured or do not come from the interaction 
region, we require leptons from J/ tp to originate from within 5 cm of the interaction point 
along the beam direction. Both lepton tracks are required to be well identified as leptons. 
Electrons are identified using a combination of specific ionization measurements [dE/dx) 
from the CDC, the ACC response, and electromagnetic shower position, shape and energy 
from the ECL Mi. Muons are identified with KLM hit positions and penetration depth Il4|. 
In order to recover dielectron events where one or both electrons have radiated a photon (final 
state radiation or bremsstrahlung) , we include the four-momentum of every photon detected 
within 0.05 radians of the original e"*" or e~ direction in the invariant mass calculation. The 
invariant mass of the candidate J/ip ^ /i^/i^(e+e~) is required to be between — 60(— 150) 
MeV/c^ and -|-36(+36) MeV/c^ of the known J/ip mass. The asymmetric mass requirements 
are due to radiative tails. The selection criteria for J/ %p are identical to those used in ref . . 

Particle identification information from the ACC, TOF and dE/dx information from the 
CDC is used to construct likelihoods Lj for each hadron type i (i = vr, K, and p). We require 
Lp/{Lp + Lk) > 0.6 and Lt^/^L^, + L^) > 0.6 to select protons and pions, respectively. For 
the prompt p candidates, tracks that are positively identified as muons or electrons are 
rejected. 

A candidates are reconstructed via the pvr" decay channel. We require that the transverse 
impact parameters of both A daughter tracks with respect to the nominal beam axis be 
greater than 0.03 cm, the z distance between the daughter tracks before constraining the 
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A vertex be less than 12 cm, and the vertex fitting be smaller than 100. The invariant 
mass of the A candidate is required to be within ±6 MeV/c^ of the nominal A mass. These 
criteria are determined to maximize Sj \/ S + where S is the number of expected signal 
events in the signal region defined below obtained from Monte Carlo (MC), and B is the 
number of expected background events obtained from sideband data. For the calculation of 
5", we assume the branching fraction of B~ J/ipKp to be 1.0 x 10~^ which is consistent 
with our final result. We apply vertex and mass constrained fits for J/ip and A candidates 
to improve the momentum resolution. 

We identify B^ J/ipKp candidates using two kinematic variables calculated in the 



center of mass frame: the beam-energy constrained mass, Mbc = y-E'beam ~'Pbi ^^^1 the 
mass difference, AM^ = Mb — itib-, where -Ebeam is the beam energy, pb and Mb are 
the reconstructed momentum and mass of the B candidate, and m^- is the nominal B~ 
mass. Candidates with \AMb\ < 0.20 GeV/c^ and Mbc > 5.20 GeV/c^ are selected for the 
final analysis. The signal region is defined to be 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ and 
|AMb| < 0.03 GeV/c^, which corresponds to three standard deviations based on the MC 
simulation. We also define a AM^ sideband region as 0.06 GeV/c^ < AM^ < 0.20 GeV/c^. 

In this analysis, we do not use the widely used AE = E's — i^beam, where Eb is the energy 
of the reconstructed B. For the J/ipAp events, the kinematic limit of AE can be expressed 



as AE = yMtot + -Eboam ~ ^hc "-^bcam, whcrc Mtot is the sum of the masses of J/tp, A and 
p. This kinematic limit is close to the signal region and depends on Mbc, which introduces a 
distortion in the Mbc — AE phase space. Consequently, events outside the Mbc/Ai? signal 
region cannot be used to model the AE/M]^^ distributions for the background. Furthermore, 
we find a large negative correlation between AE and Mbc for signal events. On the con- 
trary, AMb and Mbc are uncorrelated, which is confirmed using Monte Carlo simulations. 
Therefore, we use AM^ and Mbc to identify signal candidates in this analysis 3 ■ 

Around 15% of the selected events have more than one B candidate. We select the best 
candidate by first choosing the A candidate with the smallest total of vertex and 
mass constrained fit. If multiple B candidates with the same A remain, we select the one 
with the smallest of the J/ ipp vertex. 
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FIG. 1: (a) {AMb-, -^bc) distribution of B'^ J /il)h.p candidates, and its projections onto (b) Mbc 
and (c) AMb- The dashed lines indicates the bands used for the projections. The curves illustrate 
the result of the fit described in the text. 

Figure ^ shows a scatter plot of AMb versus Mbc and their projections for candidates 
after all selection criteria are applied. The AM5(Mbc) projection is shown for candidates 
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in the Mbc(AMB) signal region. There are six events in the (MbcAM^) signal region. 
We obtain the signal yield by fitting the Mbc distribution, since this variable has better 
agreement between data and MC than AM^. The background is described by an ARGUS 
function [13] and the signal PDF is modelled by a sum of two Gaussians plus a Crystal Ball 
line function jisjl to account for the small tail, with parameters obtained from MC. In the 
fitting, -Ebeam and the width of the main Gaussian are fixed as 5.289 GeV and 2.62 MeV/c^ 
respectively, based on a control sample of B~ D^tt~ . We simultaneously fit the signal 
and sideband regions, where the PDF in the sideband region is given by an ARGUS function 
with the same parameters as used in the signal region, except for the normalization, which is 
allowed to float. The Mbc fit gives 3.5±2.3 signal events and 0.92±0.34 background events in 
the signal box. The statistical significance of the signal, defined as y^— 2 ln(Lo/Lmax); is 2.3, 
where Lmax and Lq denote the maximum likelihood with the fitted signal yield and with the 
yield fixed at zero, respectively. We examine the contribution of other decay modes including 
J /if) and baryons, such as J/ipTfp, which may also peak in the Mbc signal region. 

Using MC, we find them to be negligible with an assumption that their branching fractions 
are comparable to B~ J/%l)A.p. As a cross check we also fit the AM^ distributions; the 
obtained signal yield is 4.7 ±2.7, with a statistical significance of 2.8, while the background 
yield is 0.55 ± 0.31. These results are consistent with the results of the Mbc fit. Since the 
signal yield is not substantial, we give an upper limit of the branching fraction as the main 
result of this analysis. 

The reconstruction efficiency (e) is estimated using signal MC where a three-body phase- 
space model is employed. We obtain e = (6.3t|3)%. The sources and amounts of systematic 
uncertainties are summarized in Table B The uncertainty of the tracking efficiency is esti- 
mated by adding the momentum dependent single track systematic error. The uncertainty 
is ~ 1% for leptons from J/ip, ~ 3.5% and ~ 1.6% for low momentum pions and protons 
from A, and ~ 1.4% for prompt protons from B. The uncertainty of the A reconstruction 
is determined by comparing the proper time distributions for data and MC simulation. For 
the uncertainty due to modeling three-body decays in the phase space, we conservatively 
assign the maximum variation of the efficiencies among the slices of M{J/ip,A), M{J/ilj,p), 
and M(A,p). 



Source 


Relative systematic error 


Tracking 
PID (proton and pion) 
Lepton identification 

A reconstruction 
Simulation modeling 
MC statistics 


±8.5% 

±8% (3% per p, 2% per vr) 
±4% (2% per lepton) 
±6% 
±41.3% / -34.1% 
±1.8% 


Total 


+43.6% / -36.8% 



TABLE I: Summary of uncertainties in the reconstruction efficiency. 

We use the method of Feldman and Cousins [lfll|. including the uncertainties of the back- 
ground and efficiency estimations [2ol|. to obtain a 90% confidence interval for the branching 
fraction given by Ns/{e x N^^ x B{J/ip —>■ x B{A pvr~)), where Ns is the signal 

yield, is the number of BB pairs, (85.0 ±0.5) x 10^, and the decay branching fractions 
B{J/ip and B{A p^r^) are taken from the world averagesfSi]. The fractions of 
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neutral and charged B mesons produced in T(45') decays are assumed to be equal. With six 
observed candidates, 0.92 ± 0.34 background events, and the uncertainties in e, N^^, and 
secondary branching fractions mentioned above, we obtain an upper bound of the interval 
from this procedure of B{B^ J/ipAp) < 4.1 x 10~^, which we interpret as a conservative 



estimate of the 90% confidence upper limit of the branching fractional 
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FIG. 2: Momentum distributions of (a) J/ip, (b) A and (c) p in the rest frame of the reconstructed 
B for the six B~ — > J /ij)Kp candidates. The histograms are phase space distributions from signal 
MC normalized to six events. 

Figure 121 shows the momentum distributions for J/ip, A, and prompt p of the six candi- 
dates in the B rest frame. We do not observe any significant enhancement above the phase 
space distribution, such as would be expected by an intermediate resonancef^. 

In summary, we have searched for the decay of B~ — > J/ipKp at Belle with 78 fb~^ data 
collected at the T(4S') resonance. No statistically significant signals are found. We set 
an upper limit of B{B~' J/ipAf)) < 4.1 x 10^^ at 90% confidence level. This result is 
consistent with the BaBar result i9]. This mode does not account for a significant fraction 
of the observed excess in the low momentum region of 5 — > J/ipX. 
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